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Rapidly spinning massive black holes in active galactic nuclei: 
evidence from the black hole mass function 
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ABSTRACT 

The comparison of the black hole mass function (BHMF) of active galactic nuclei (AGN) 
relics with the measured mass function of the massive black holes in galaxies provides strong 
evidence for the growth of massive black holes being dominated by mass accretion. We derive 
the Eddington ratio distributions as functions of black hole mass and redshift from a large 
AGN sample with measured Eddington ratios given by KoUmeier et al. We find that, even at 
the low mass end, most black holes are accreting at Eddington ratio A ^ 0.2, which implies 
that the objects accreting at extremely high rates should be rare or such phases are very short. 
Using the derived Eddington ratios, we explore the cosmological evolution of massive black 
holes with an AGN bolometric luminosity function (LF). It is found that the resulted BHMF 
of AGN relics is unable to match the measured local BHMF of galaxies for any value of 
(constant) radiative efficiency ryrad- Motivated by Volonteri, Sikora & Lasota's study on the 
spin evolution of massive black holes, we assume the radiative efficiency to be dependent 
on black hole mass, i.e., 77iad is low for A/bh < IO^Mq and it increases with black hole 
mass for Mbh ^ IO^Mq. We find that the BHMF of AGN relics can roughly reproduce the 
local BHMF of galaxies if Ty^ad ^ 0.08 for Mbh < lO^M© and it increases to > 0.18 for 
Afbh ^ IO^Mq, which implies that most massive black holes (> IO^Mq) are spinning very 
rapidly. 

Key words: (galaxies:) quasars: general — accretion, accretion discs- 
galaxies: evolution 
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1 INTRODUCTION 

It is believed that quasars are powered by accretion on to massive 
black holes, and the growth of the massive black holes could be 
governed by mass accretion in quasars. The massive black holes 
(AGN relics) should be present in the centres of galaxies. Thus, 
the luminosity functions (LF) of active galactic nuclei (AGN) pro- 
vide important clues on the growth of massive black holes. It 
is indeed found that most nearby galaxies contain massive black 
holes at their centres, and a tight correlation is revealed between 
central massive black hole mass and the velocity dispersion of 
the galaxy ( .Ferrarese & Me rritt 2000; Gebhardt_et al, 2000l) . The 
black hole mass is also found to be tightly correlated with the 
luminosity of the spher oid component of it s host galaxy (e.g., 
iMagorrian et al . 1998; M arconi & Hunlll2003l) . These correlations 
of the black hole mass with velocity dispersion/host galaxy lumi- 
nosity were used to derive the mass functions of the central massive 
black holes in galaxies (e.g., Yu & Tremaind 20021 ; Marconi et al.l 



can also be calculated by integrating the continuity equation of 
massive black hole number density on the assumption of the growth 
of massive black holes being dominated by mass accretion, in 



of AGN (e.e.. ICavaliere. Morrison & Wood 197 ll: SoltaiJ 19821; 


Chokshi & Turner 19921; 


Small & Blandfordj 1991 Marconi et al.1 


2004;iShankaretalJl2004 


;lTamura. Ohta & Ueddbooel). Such cal- 



culations on the cosmological evolution of massive black holes 
were usually carried out by adopting two free parameters: the ra- 
diative efficiency 77i.ad and the Eddington ratio A for AGN. The 
derived BHMF of AGN relics in this way is required to match 
that estimated either from velocity dispersion or the luminosity 
of the spheroid of its host galaxy, which always leads to A ~ 1, 
almost all AGN are required to be accreting close to the 



Eddington liinit (eg.. TYu & Tremaind 1 20021: [Marconi et ani2004 



|2004 iTamura. Ohta &Uedal Hood lGraham etalJ l2007h . On the 
other hand, the black hole mass function (BHMF) of AGN relics 
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IShankar et al.ll2004l ; Tmnura. Ohta & Uedj|2006h 

In the last decade, several approaches for measuring the 
masses of the central black holes in AGN have been developed, 
in which the reverb eration mapping m ay be the most effective 
one (Petersorj |T993l ;[ Kaspi et al. 2000). Using the tight correla- 
tion between the size of the broad-line region and the optical lu- 
minosity established with the reverberation mapping method for 
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a sample of AGN, the black hole masses of AGN can be easily 
estimated from their optical luminosity and width of broad emis- 
sion line. The Eddington ratios for thousands of AGN were esti- 
ma ted with the analys e s of th e Sloan Digital Sky Survey (SDSS) 
bv iMcLure & Dunlod J2004h . which indicate that the mean Ed- 
dington ratio -Lboi/J/Edd — . 1 at z ~ 0.2 to ~ 0.4 at z ~ 2. 
IWamer. Hamann & Dietrichetl ( |2004|) also derived the Eddington 
ratio distribution for a sample of ^ 50 AGN with redshifts < 
z < 5. As pointed bv ^ Kollmeier et all ( l2006h . both these derived 
Eddington ratios are heavily weighted towa rds high-luminosity ob- 
jects due to the limited sensitivity of SDSS. lKoUmeier et alj jiool) 
estimated the Eddington ratios of AGN discovered in the AGN and 
Galaxy Evolution Survey (AG ES), which is more sensitive than the 
SDSS jKochanek et alj|2004h . The derived Eddington ratio distri- 
bution at fixed luminosity is well described by a single lognormal 
distribution pe aked at ^ 0.25 independent of redshift and luminos- 
ity (see .KoUmeier et al.ll2006l . for the details). 

In this w ork, the Eddington rati o distribution at fixed lumi- 
nosity given bv'KoUmeie r et al. is converted to that at fixed 
black hole mass by using an AGN LF. We integrate the continuity 
equation for black hole number density adopting the derived Ed- 
dington ratio distributions of AGN to calculate the BHMF of AGN 
relics at different redshifts z, which is different from a free pa- 
rameter A adopted in most previous works. The resultant BHMF 
of A GN relics is constrained by those estimated from the galaxy 
LFs jMarconi et al. I l2004l : iTamura. Ohta & Uedj|2006h . The con- 
ventional cosmological parameters Qm = 0.3, Qa = 0.7, and 
Ho = 70 km Mpc~^ have been adopted in this work. 



2 THE EDDINGTION RATIO DISTRIBUTION OF 
ACTIVE GALACTIC NUCLEI 

The Eddington ratio distribution of AGN for given bolometric lu- 
minosity can be approximated as a log-normal distribution: 



2Tva 



exp 



ii-i^r 



2a^ 



(1) 



where l = log A, A = iboi/iEdd, M — log 0.25 and cr ~ 0.3 (see 
iKollmeier et ai]|2006l. for the details). We can derive the BHMF of 
AGN from the bolometric LF $(2, Lboi): 



A*'agn(z, Mbh) 



<I>(z,Lb 



^ d log Lboi 
d log Afbh 



(2) 



where log Lboi = I + logLEdd, ^Edd = 1.251 x 
lO'^^Mbh ergs~^, and Mbh is in units of solar mass. Using the 
bolometric LF <l?(z, Lboi) of AGN, the Eddington ratio distribution 
for given black hole mass Mbh can be calculated with 



x(z,Mbh,0 



$(z,Lboi)C(0 

iVAGN(z,Mbh)' 



(3) 



where Lboi = lO'LEdd, and the BHMF of AGN, Nagn{z, Mbh), 
is available with Equation lO. The mean Eddington ratio of AGN 
with Mbh at z is 



A(2,Mbh)= / Ax(z,Mbh,Od^- 



(4) 



In this work, we adopt the luminosity-dependent density evo- 
lution (LDDE) bolometric LF calculated from the rest-frame opti- 
cal, soft and hard X-r ay, and near- and mid-IR bands in the re dshift 
interval 2 = — 6 bv lHopkins, Richards & Hemguistl ( |2007|) : 

$(Lboi, z) = $(Lboi, 0)ed(Lboi, z) 




Figure 1. The Eddington ratio distributions for given black hole mass Afbh- 
The red, green, and blue lines represent z = 0, 1, and 3, respectively. The 
solid, dashed, and blue lines are for Mbh = , 10^, and 10^ Mq, re- 
spectively. 



(Lbol/L*)Tl + (Lbol/L,)T2 

The density function ed is given by 



ed(Lboi, z) 



ed(Lboi, z) = 
with 

Zc(Lbol) 

and 

pl(LboO 



pi 



(5) 



Z ^ Zc 



(1 + Z,r'[{l + Z)/(1 + Zo(Lbol))]^' Z > Zo 



Zc,o(Lbol/Lc)" (Lboi ^; Lc), 
Zc,0 (Lboi > Lc), 



pl46 + Pi [log(Lboi/10'"'ergs s ^) 
p2(Lboi) =p2i6 + /32[log(Lboi/10''^ergs s"^) 



(7) 

(8) 
(9) 



All the parameters of the LF are as follows: log (f)* = —6.20 ± 
0.15 Mpc"^ log L. (ergs s-i) = 45.99 ± 0.10, 71 = 0.933 ± 
0.045, 72 = 2.20±0.14, log Lc(ergs s'^) = 46.72±0.05, Zc,o = 
1.852 ± 0.025, a = 0.274 ± 0.025, pUe = 5.95 ± 0.23, _p246 = 
-1.65 ± 0.21, /3i = 0.29 ± 0.34, and = - 0.62 ± 0.17 (see 
Table 4 in lHopkins. Richards & Hemquistll2007h . 

In Fig. [1] we plot the Eddington ratio distributions 
x(z, Mbh, I) for fix ed black hole mass deri ved from that for fixed 
luminosity given by iKoUmeier et alj ( I2OO6I) . We find that the de- 
rived Eddington ratio distributions are close to the lognormal dis- 
tribution, while their peaks vary with black hole mass and redshift. 
We plot the mean Eddington ratios as functions of black hole mass 
in Fig.|2] It is found that the mean Eddington ratios A(z, Mbh) are 
in the range of ~ 0.1 — 0.3 as functions of black hole mass Mbh 
and redshift z. 



3 THE EVOLUTION OF MASSIVE BLACK HOLES 

Th e evolution of massive bla ck hole number density is described 
bv dSmall & Blandforj|l992h 
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The AGN LF plays an important role in the study of the cos- 
mological evol ution of massive black holes. The optical quasar LF 
was ado pted inlYu & Tremainel l l2002h . however, the optical quasar 
LF (e.g.. lBovleet altood) has missed faint AGN (either intrinsic 
low-luminosity or obscured AGN). The hard X-ray surveys (~ 2 — 
lOkeV) can trace the whole AGN population, including obscured 
type II AGN. The hard X-ray LF derived by Ueda et al. qoolh 
was used in some works on the cosmological evolution of mas- 



10 



dt 



+ 



sive black holes (e.g., 'Ma rconi et al.ll2004l : lTamura. Ohta & Uedal 
bOOdlShankar et al. 2007). The number density of Compton-thick 
AGN is still quite uncertain, which is not included in the hard X- 
ray LF. The contribution of Compton-thick AGN to the black hole 
evolution was taken into account by multiplyi ng a correction factor 
of 1.6 independent of the luminosity (e.g., see lMarconi et alj2004l : 
iTamu ra. Ohta & Ueda 2006). 

The hard X-ray (> 20 keV) and the mid-IR (5 - 50 /im) 
bands are optim al for detection of AG N with column densities < 
10 cm" (e.g., lTreister & Urrvl2005l) . The observations with the 
International Gamma-Ray Astrophysics Laboratory (INTEGRAL) 
and the Swift Burst Alert Telescope (BAT) indicate that the frac- 
tion of th e absorbed AGN decreases with the 20 — 100 keV lu- 
minosity iMarkwardt et alJl2005l : lB^assani et ^l2006l) (but also see 
IWang & Jiang 2006h . which is confirmed by the mid-IR Spitzer 
observation s of 25 lu minous and distant quasars ( Maioli no et al.l 
|2007). Ma iolino et al.l ( 2007) suggested that the fraction of the ob- 
[iV(Mbh, t) < M(Mbh, t) >] = ^(Afbh, i), (10) scured AGN to the total can be well fitted with 



10° 



10^ 



10" 



Figure 2. The mean Eddington ratios as functions of black hole mass M^h 
for the AGN at different redshifts: z = O(black), l(red), 2(green), 3(blue) 
and 4(yellow). 



9iV(Mbh,t) d 



where N{Mhh, i) is the mass function of massive black holes in- 
cluding both active and inactive black holes, < M(Mbh,t) > is 
the mean mass accretion rate for the black holes with Mbh, and 
^(Mbh, i) describes the effect of black hole mergers on the BHMF. 
The total black hole mass density will not be altered by merg- 
ers, if t he mass loss caused by the gravitational radiation is ne- 
glected. ^hankar et al J ( 1200 7*) assessed the importance of black hole 
mergers on the evolution of the BHMF using a simple mathemat- 
ical model that assumes constant probability Pmer g of equal mass 
merge rs per Hubble time, similar to the models of Richstone et al.l 
( Il998l) . They found that the effect of black hole mergers is to 
slightly lower the number density of small black holes and in- 
crease the number density of massive black ho les, if Pmerg = 0.5 
is adopted (see Fig. 13 in IShankar et al]|2007l) . Observational es- 
timates of the galaxy merger rate and its mass dependence span a 
substantial rang e (e.g.. Bell et al. 2o63; IConrov. Ho & Whitel2007l : 
iMas iedi. Hogg & Blanton 2008), and Pmcrg = 0.5 is roughly con- 
sistent with the high end of these estimates (see lShankar et alj2007l 
for the detailed discussion). Thus, they concluded that the impact 
of mergers on the BHMF is rather small compared with mass accre- 
tion. In this work, we neglect the effect of black hole mergers, i.e., 
'S'fMbh, t) = in Eg . {TOt. as in most of the previous works (e.g.. 



/o 



where 



•^ODt 



1 4_ r0.414 ' 



(14) 



opt 



ALa(5100A)[: 



erg s 



. fl.J . 

Yu & Tremai"n3l2002l : [Marconi et aLlliooi ; iTamura. Ohta & Ued3 
2006I) . The mean mass accretion rate is 



< M(Mbh,t) >= ^(Ml±l^)^5(Mbb,t), 
where the duty cycle of active black holes is defined as 

iVAGN(Mbh,t) 



(11) 



(12) 



Substituting Eqs. Jilt and il2\ into Eq. dlOb . we can rewrite the 
black hole evolution equation as 



]^q45.63 

iMtiller & Hasinged ( |2007|) found that the fraction of type II 
AGN detected in the hard X-ray band can be described by this 
function (Equation II4t quite well (s ee Fig. 3 in their paper). 
iHopki ns. Richards & Hemguistl ( l2007h suggested that the fraction 
of Compton-thick to the total also decreases with luminosity and 
it is less than ~ 30 per cent based on a variety of very hard 
X-ray/soft gamma-ray observations on AGN (see their paper for 
the detailed discussion). Motivated by the results of these works, 
besides the luminosity-independent correction for Compton-thick 
AGN, we tentatively emp loy a similar l u minos ity-dependent cor- 
rection as that given by .Maiolino et alj ( l2007h . We assume the 
number ratio of Compton-thick to Compton-thin AGN to be 

/CT = . , °^0.414 ' (15) 
^ '-hoi 

wher e £boi = Lbol/10 ^^''■^ as Lboi ==; 9ALa(5100A) is adopted 
(e.g., lKasDiet"Ii]|2000l) . We change the numerator in Eq. il4\ to 
0.6 here, so that /ct redu ces to ^ 0.6 for low- luminosity AGN, 
which is the same as that in lMarconi et al. 

|i2004), while /cT 

for luminous AGN. 

In most of the previous works, both the radiative efficiency of 
r/raci and Eddington ratio A are free parameters, and the compar- 
isons between the BHMF of AGN relics and the measured local 
BHMF of galaxies always require: T^rad ~ 0.1 and A ~ 1 (e.g., 



Yu & Tremain3 12002| : iMarconi et alj|2004 ITamura. Ohta & Uedj 
20061) . As we have derived the mean Eddington ratio distributions 



dNjz, Mbh) 
dz 



dt d 



A(2, Mhh)LEddNAGN{z, Mb 



.(13) 



as functions of black hole mass and redshift in the last section, there 
is only one free parameter r/^ad in our calculations for the cosmo- 
logical evolution of massive black holes. The local BHMFs esti- 
mated by using the correlation of the black hole mass with host 
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galaxy l uminosity are adopted in this work (see, iMarconi et alj 
[2004 ; Tamura. Ohta & Uedd[20o3 . for the details). The continu- 
ity equation l |13l l for black hole number density is integrated from 
z = «max by using Eqs. l|2j-l|4j and assuming the duty cycle is 
0.5 at Zmax- The final results are insensitive to the initial condi- 
tions at Zmax. In all our calculations, Zmax = 4 is adopted, because 
the Eddington rati o distributions are cal culated from a sample of 
AGN with z < 4 jKollmeier et"aLll2006h . The resulted BHMFs of 
AGN relics at low redshifts are insensitive to the value of 2;max, 
because the fraction of local black hole mass accreted at high red- 
shifts can be neglected. We plot our results with different val- 
ues of r^rad in the upper panel of Fig. |3] which indicates that the 
measured local BHMF cannot be fitted with any values of rji^d- 
This is due to the mean Eddington ratios derived in this work be- 
ing 0.1 — 0.3, which deviates significantl y from A ~ 1 sug- 
gested in most of th e previous works (e.g.. IYu & Trem^n3l2002l : 
iMarconi et al ] |2004h . In our calculations, we only consider the un- 
certainty of the number density in th e bolometric LF given by 
iHopkins, Richards & Hemguisll fioOTh . 

The radiative efficiency of black hole accretion is closely re- 
lated to the black hole spin. The massive black holes will be spun 
up through accretion, as the black holes acquire mass and angu- 
lar momentum simultaneously through accretion. The spins of the 
massive black holes may also be affected by mergers of black holes. 
A rapidly rotating new black hole will be present after the merger of 
two black holes, only if the binary's larger member already spins 
quickly and the merger wit h the smaller hole, or if the binary's 
mass ratio approaches unity jHughes & Blandforj|2003h . The co- 
moving space density for heavier black holes is m uch lower than 
that f or smaller black holes (e.g., see the BHMF in IMarconi et al.l 
l2004h . which means that the probability of the mergers of two 
black holes with similar masses is lower for heavy black holes. 
This implies that the spins of heavy black holes are mainly regu- 
lated by accretion rather than the mergers, and the spin parameters 
can reach ~ 1 after their masses are doubled through accretion. 
For smaller black holes in disk galaxies, a small number of minor 
mergers might have happened, which are believed to be responsible 
for rebuilding their host galaxy disks, and small minor accretion 
episo des on to black holes are trigger e d by these minor mergers 
(e.g., l\^lonteri. Sikora & Lasotal uOOtI . iKing. Pringle & HofmannI 
I2OO8I) . This is also supported by the observations that single ac- 
cretion events last ~ 10^ years in Seyfert galaxies and their 



total active lifetime is 
Ho. F ilippenko & SargenlJ_1997 
Volonteri. Sikora & Lasot; 



; 10" 



10" 



years 



CKharb et alj l2006l: 



ot32 



. .Volonteri. Sikora & Lasotal2007r) . 
( l2007h~ studied on how the accretion 
from a warped disc influences the evolution of black hole spins 
and concluded that within the cosmological framework, one in- 
deed expects most supermassive black holes in elliptical galaxies 
to have on average higher spin than black holes in spiral galax- 
ies. The random small accretion episodes (e.g., tidally disrupted 
stars, accretion of molecular clouds) might have played a more im- 
portant role on the spin evolution of small black holes in spiral 
galaxi es, which lead to relatively low average spins for these black 
holes jVolonteri. Sikora & Lasot j|2007l : [King. Pringle & HofmannI 
I2OO8 I). Motivated by their results on the cosmological spin evolu- 
tion of black holes, we tentatively adopt a i\/bii-dependent radiative 
efficiency, in which r/rad remains constant for A/bh < 10* and in- 
creases as a power-law with black hole mass for Afbh J? 10* : 




r/rad 



if Afbh < IO^Mq; 
if A'/bh ^ 1O*M0 



(16) 



log M^, 

Figure 3. The local black hole m ass functions. The gr een dotted line repre- 
sents the local BHMF derived by Mai'coni et al.' 1 2004) (see their paper for 
the details). The shaded area is the BHMF tr ansformed from the COMBO - 
17 LF by using Mbh - L relation (taken from lTamura, Ohta & Uedal20O6l) . 
The BHMFs of AGN relics calculated from the AGN LF are plotted in 
the upper panel for different radiative efficiencies (ry^ad = 0.05, 0.1, and 
0.2, from up to down). In the lower panel, we calculate the BHMFs of 
AGN relics from the AGN LF using A/bh-dependent radiation efficiency 
^rad(Afbh) (see the text for the details). The blue lines are calculated with 
a correction factor of 1 .6 in dependent of luminosity for the Compton-thick 
AGN, while the red lines represent the cases with a luminosity-dependent 
correction for the Compton-thick AGN (see the text for the details). 



in our calculations. We find that the measured local BHMF can be 
roughly reproduced by the BHMF of AGN relics provided J7rad,o ~ 
0.08 and g = 0.36 are adopted (se e the lower panel in Fig.[3]l. 

iTamura. Ohta & Uedal ( l2006h tried to derive the BHMFs with 
redashifts up to z ~ 1 from the spheroid LF of early-type galaxies 
using the correlation between the spheroid luminosity and black 
hole mass (see their paper for the details), which provide further 
constraints on the model calculations for the cosmological evolu- 
tion of massive black holes. The model calculations performed with 
this A/bh -dependent radiative efficiency (Equation |16t are com- 
pared with the spheroid-BHMFs derived by Tamura, Ohta & Ue(3 
120061) in Fig. g] We find that calculated BHMFs of AGN relics 
can roughly reproduce the spheroid-BHMFs at different redshifts 
either for luminosity-dependent or luminosity-independent correc- 
tions for the Compton-thick AGN. 



4 DISCUSSION 

As in the most previous works, we implicitly assume that the black 
hole growth is dominated by mass accretion in bright AGN, while 
some inactive black holes may still be accreting gases, though their 
mass accretion rates are very low. If the duration of the accretion 
in these objects is as long as the Hubble timescale, they can ac- 
crete sufficient mass comparable with that accumulated in bright 
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Figure 4. The BHMFs at different redshifts of 2 = 0, 0.25, 0.45, 0.65, 0.85, 
and 1.05 transfo mied from th e COMBO-17 LFs are indicated with shaded 
ai-eas (taken from lTamuraTohta & Ueda 2006) . The BHMFs of AGN rehcs 
at different redshifts calculated from the AGN LF are plotted by using Af^h" 
dependent radiation efficiency r)^s.d (Afbh) (see the text for the details). The 
blue lines are calculated with a correction factor of 1.6 in dependent of 
luminosity for Compton-thick AGN, while the red fines represent the cases 
with a luminosity-dependent correction for the Compton-thick AGN. 



AGN phases, as the AGN phase is much shorter than the Hub- 
ble timescale. It is believed that the advection dominated accre- 
tion flows (ADAFs) are present in those objec ts, which are ver y 
hot and radiate mostly in hard X-ray bands (Naray an & Yilll994l) . 
They are very difficult to be detec ted due to low luminosity, un- 
less those in the nearby Universe. ICad jZOOSh suggested that the 
accretion of such low-luminosity objects can be constrained by the 
hard X-ray background, though the emission from most of these 
individuals cannot be detected by any facilities now. It was found 
that less than ~ 5 per cent of the local black hole mass density 
was accreted during the ADAF phases, whi ch will be even lower 
if the Compton-thick AGN are included (seelCaoll2007l for the de- 
tails). |Ho2yns[^^ayan|&^emguis3 1 I2OO6I) considered the distri- 
bution of local supermassive black hole Eddington ratios and ac- 
cretion rates, accounting for the dependence of radiative efficiency 
and bolometric corrections on the accretion rate. They also found 
that black hole mass growth was dominated by AGN phase, and 
not by the radiatively inefficient low accretion rate phase in which 
most local supermassive black holes are currently observed. 

The main difference of this work from the previous works is 
that the Eddington ratio distributions are derived from an AGN 
sample with measured Eddington ratios jKolImeier eTall l2006h . 
The Eddington ratio distributions for fixed black hole mass derived 
in our work approximate to the lognormal distribution (see Fig.[T]l, 
and the mean Eddington ratios are in the range of ~ 0.1 — 0.3 
varying with black hole mass and redshift (see Fig. |2j. For most 
cases, the mean Eddington ratios peak at ~ IO^Mq, and then 
decline with increasing black hole mass. Even at the low mass 
end, most black holes are accreting at A ~ 0.2, which implies 



that the objects accreting at extremely high rates should be rare 
or such phases are very short. It was suggested that the radiative 
efficiency r^rad declines for a slim accretion disc provided the mass 
accretion rate is sufficiently hig h due to the photon trapping effor t 
(e.g.,[A bramowicz et al. 1988; B egelman|[T973 : IWang et al.lfl999h . 
IWatarai et al.l (12000)'s calculations on the slim discs showed that 
the radiative efficiency will not deviate significantly from that for 
standard thin discs if Lboi/iEdd ^ 2, which implies that the 
present adopted radiative efficiency independent of Eddington ratio 
A is indeed a good assumption. 

There is only one free parameter ry^ad in our calculations for 
the cosmological evolution of massive black holes. We find that 
the resulted BHMF of AGN relics is unable to reproduce the mea- 
sured local BHMF for any value of rj,.a_d adopted, provided the ra- 
diative efficiency r^rad i s independent of black hole mass, as treated 
i n previous works (e.g., | Yu & Tremainell2002l : iMarconi et alj2004 ' 



iTamura, Ohta & Uedall200^ The mean Eddington ratios adopted 
in our calculations are derived from an AGN sample, which are 
in the range of ~ 0.1 — 0.3. Thus, it is not surprising that the 
local BHMFs cannot be reproduced by our calculations with any 
constant radiative efficiency, because the Eddington ratio A ~ 1 
is usually required in order to let the resulted BHMF match the 
local one in those works. In this work, we use two different cor- 
rections (either luminosity-independent or luminosity-dependent) 
for the Compton-thick AGN (see Sect. 3 for the details), and find 
that the final results are quite similar (see Figs.[3]and|4l. We also 
use the hard X-ray LF derived from an AGN sample at high red- 
sh ifts by [Si lverman et al. (200^ i n stea d of the bolometric LF 
of iHopkins. Richards & Hernguis^ in the calculations. It is 

found that the main results of this work change very little and the 

mai n conclusion is not altered. 

IVolonteri. Sikora &. Lasotal ^20071) studied on how the accre- 
tion from a warped disc influences the evolution of black hole spins 
and concluded that within the cosmological framework, one indeed 
expects most supermassive black holes in elliptical galaxies to have 
on average higher spin than black holes in spiral galaxies, where 
random, small accretion episodes (e.g., tidally disrupted stars, ac- 
cretion of molecular clouds) might have played a more important 
role. Thus, we tentatively adopt a A/bh-dependent radiative effi- 
ciency (see Eq.|I6l>, in which r/^ad remains constant for Mbh ^ 10** 
and increases with black hole mass for A/bii > 10*. This Mbh- 
dependent radiative effic iency is qua li tative ly consistent with the 
results of Volo nteri. Sikora & Lasotal l l2007h . It is found that the 
measured BHMFs can be fairly well reproduced by our model cal- 
culations with this Mbh-dependent radiative efficiency (see Figs.|3] 
andgll, which require ry^ad > 0.18 for Mbh > IO^Mq. This pro- 
vides evidence for most massive black holes being spinning very 
rapidly. It is interesting to find that a ~ 0.9 is also required by 
the fitting of the residual hard X-ray background wi th the emis - 
sion from the ADAFs in the low-luminosity objects l lCaoll2007h . 
Our calculations can be improved if the mean spin parameter a as 
a function of black hole mass is available from the work within 
the cosmological framework, which is beyond the scope of present 
work. 

In our present calculations of the black hole evolution, the 
black hole mergers have been neglected. Shankar et al. ( 2007) as- 
sessed the importance of the black hole mergers on the evolu- 
tion of the BHMF. They found that the impact of black hole 
mergers on the cosmological evolution of BHMF may prob- 
ably be small compared with black hole accretion processes, 
while its im pact on the black hole spin evolution may be im- 
portant (e.g.. IWilson & Colbertlll995l ; iHughes & BIandfordll2003l: 
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IVolonteri et al ] |2005l : IVolonteri, Sikora & Lasota"2007') . The effect 
of black hole me rgers increases the n umber density of very mas- 
sive black holes iShankar et al.ll2007h . which implies that the ra- 
diative efficiencies for very massive active black holes should be 
higher than the present values if black hole mergers are included in 
our calculations. This strengthens our conclusion that most massive 
black holes are spinning very rapidly. 
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